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Abstract 

We consider a sterile neutrino to be an unparticle, namely an unsterile neutrino^ with anomalous 

dimension 77 and study its mixing with a canonical active neutrino via a see-saw mass matrix. We 

show that there is no unitary transformation that diagonalizes the mixed propagator and a field 

redefinition is required. The propagating or "mass" states correspond to an unsterile-like and 

active-like mode. The unsterile mode features a complex pole or resonance for < 77 < 1/3 with an 

"invisible width" which is the result of the decay of the unsterile mode into the active mode and 

the massless particles of the hidden conformal sector. For 77 > 1/3, the complex pole disappears, 

merging with the unparticle threshold. The active mode is described by a stable pole, but "inherits" 

a non-vanishing spectral density above the unparticle threshold as a consequence of the mixing. 

We find that the radiative decay width of the unsterile neutrino into the active neutrino (and a 

r 1 ^ 

photon) via charged current loops, is suppressed by a factor ^ 2sin^(0o) ^ ^ where is the 
mixing angle for 77 = 0, M is approximately the mass of the unsterile neutrino and A ^ M is 
the unparticle-scale. The suppression of the radiative (visible) decay width of the sterile neutrino 
weakens the bound on the mass and mixing angle from the X-ray or soft gamma-ray background. 

PACS numbers: 98.80.Cq;12.60.-i;14.80.-j 



* Electronic address: "boyan@pitt.edu 

"f" Electronic address: rh4a@andrew.cmu.edu 

■'■Electronic address: j hut asoi@ andr e w . emu . edu 



1 



I. INTRODUCTION 



Neutrino masses and oscillations are now an established phenomenon and an undisputable 
evidence of physics beyond the standard model. Although the origin and scale of masses 
remains a challenging question, the see-saw mechanism provides a compelling explanation 
of small active neutrino masses as the result of ratios of widely different scales [Ij. Among 
extensions of the standard model, the addition of sterile neutrinos, namely SU(2) singlets 
with a mass in the few keV range, acquires particular importance as a potential warm dark 
matter candidate [SElllElEllZllHlinilinillllliailS could provide possible solutions 

to a host of astrophysical problems [6]. The radiative decay [15] of a sterile-like neutrino 
mass eigenstate into an active-like mass eigenstate and a photon leads to a decay line that 
could be observable in the X-ray or soft gamma ray background [16j. The non-observation 
of this line provides a constraint on the mass and mixing angle of sterile-like neutrinos 

More recently sterile neutrinos with mass in the GeV range have been proposed as 
explanations of two seemingly unrelated and unexpected phenomena: an excess of air 
shower events at the SHALON gamma ray telescope, in a configuration where the ex- 
pected number of events is negligible pOj, and as a potential explanation of the MiniBoone 
anomaly [21j, namely the prominent peak of electron-neutrino events above background for 
300 MeV < Ey < 475 MeV. Sterile neutrinos in a three-active, two-sterile (3 + 2) oscillation 
scheme were proposed in Ref. [22j as a possible explanation of the MiniBooNE anomaly, 
and an alternative explanation invoking the radiative decay of a heavy sterile neutrino with 
a small magnetic moment was proposed in Ref. [23] . 

In this article we study the possibility that sterile neutrinos are a manifestation of un- 
partides^ which then mix with the active neutrinos via a see-saw-type mass matrix. 

In a recent series of articles, Georgi [24j suggested an extension of the standard model in 
which particles couple to a conformal sector with a non-trivial infrared fixed point acquiring 
(large or non-pert urbative) anomalous dimensions with potentially relevant consequences, 
some of which may be tested at the Large Hadron Collider (LHC) |25l [26l [2?J. Early 
work by Banks and Zaks provides a realization of a conformal sector emerging from 
a renormalization fiow toward the infrared below an energy scale A through dimensional 
transmutation, and supersymmetric QCD may play a similar role [29j. Below this scale 
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there emerges an effective interpolating field, the unparticle field, that features an anomalous 
scaling dimension [24j. 

Various studies recognized important phenomenological [24l [25l [30l IHU [32] astrophysical 
[331 [SI ESj and cosmological [36l[3Zl[38l[39l[40l|4ll[4^ consequences of unparticles, 
including Hawking radiation into unparticles [44], aspects of CP- violation [45j, flavor physics 
[46] and low energy parity violation [47j. More recently, the consequences of mixing of 
unparticle scalar fields and a Higgs field were studied in Ref. [48j with implications for 
slow-roll inflation. 

A deconstruction program describes the unparticle as a tower of a continuum of exci- 
tations ^J. However, this is not the only interpretation of unparticles. As mentioned in 
Ref. [48], anomalous scaling dimensions are ubiquitous in critical phenomena near an in- 
frared fixed point [50l [5TJ , which is the main observation in Ref. [28j . There is also a well 
known phenomenon in QCD [52j where anomalous dimensions emerge from the multiple 
emission and absorption of gluons as a result from the resummation of infrared Sudakov 
logarithms. Similarly, in QED, Bloch-Nordsieck resummation of infrared divergences arising 
from multiple emissions and absorptions of photons yield threshold infrared divergences and 
lead to a renormalized electron propagator that also features anomalous dimensions [53l [54] • 
The renormalization group resummation of absorption and emission of massless quanta lead 
quite generally to anomalous dimensions in the propagators [55j . 

A recent study [56] suggests a connection between unparticles and the Miniboone 
anomaly, namely, that the heaviest mass eigenstates corresponding to the mixed — 
decays into the lightest eigenstate and a scalar unparticle (see also [57j). Furthermore, un- 
particle contributions to the neutrino-nucleon cross section and its influence on the neutrino 
flux expected in a neutrino telescope such as IceCube have been reported in Ref. [58j . 

If heavy sterile neutrinos decaying into lighter active neutrinos and massless particles of 
a hidden conformal sector described by scalar unparticles are a plausible explanation of the 
MiniBooNE anomaly, then the coupling of the sterile neutrinos to this degree of freedom 
will necessarily lead to the consideration of the sterile neutrino itselj being an unparticle. 

This consideration emerges naturally from the "deconstruction" argument [49j , since the 
coupling of the sterile neutrino to the scalar unparticle leads to a spectral representation of 
the sterile neutrino propagator that features anomalous scaling dimensions. Alternatively, 
(but equivalently) the emission and absorption of massless (conformal) quanta lead to in- 
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frared threshold divergences akin to Sudakov logarithms whose renormalization group (or 
Bloch-Nordsieck [53]) resummation leads to anomalous dimensions j 52l [55] . 

Massive fermionic unparticles with a soft conformal breaking mass term had been in- 
troduced in Ref. [59j and a very interesting proposal in which right handed neutrinos are 
fermionic unparticles was considered in Ref. [60j . 

If unparticle physics proves to be an experimentally relevant extension of the standard 
model, it is natural to consider sterile neutrinos, namely SU{2) singlets, as being unparticles, 
and we refer to them as unsterile neutrinos. 

Unsterile neutrinos are assumed to couple to a "hidden" conformal sector beyond the 
standard model and acquire a (possibly large) anomalous scaling dimension below a scale 
A of dimensional transmutation at which the infrared fixed point of the conformal sector 
dominates the low energy dynamics. 

In this article we consider such a possibility and study the consequences of unsterile 
neutrinos mixing with active neutrinos via a typical see-saw mass matrix. We consider the 
simplest scenario of one unsterile and one active Dirac neutrino to establish the general 
consequences of their mixing. Our objectives in this article are two-fold: 

• Because unsterile neutrinos feature non-canonical kinetic terms, novel aspects of mix- 
ing phenomena emerge. We explore the fundamental aspects of mixing between these 
unsterile and the usual active neutrinos via a see-saw type mass matrix. 

• We also focus on potential cosmological consequences, in particular if and how the 
unparticle nature of a sterile neutrino modifies its radiative decay into an active-like 
neutrino. This decay rate is an important ingredient to establish bounds on masses 
and mixing angles from the cosmological X-ray or soft-gamma background in the case 
when the mass of the sterile-like neutrino is in the keV range, which is of interest when 
considering it as a dark matter candidate. 

Our results can be summarized as the following: 

• The spinor nature of the unparticle field introduces novel aspects of mixing, in which 
there is no unitary transformation purely in fiavor space that diagonalizes the full 
propagator. The diagonalization requires a non-unitary transformation and field re- 
definition followed by a momentum-dependent transformation that is unitary below 
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the unparticle threshold but non-unitary above. The resulting mixing angles depend 
on the four-momentum^. 

• For a see-saw type mass matrix that mixes the unsterile and active neutrino, we 
find an unsterile-like mode with an "invisible" decay width. A renormalization group 
inspired "resummation" argument suggests that this width is a result of the decay of 
the unsterile-like mode into the active-like mode and particles in the "hidden sector" . 
A complex pole for the unsterile mode exists only for < < 1/3, where r] is the 
unparticle anomalous dimension. Asrj ^ l/3~ the real part of the pole approaches the 
unparticle threshold from above and the width becomes large. For r] > 1/3 the spectral 
density for the unsterile mode does not feature a complex pole, but is described by a 
broad continuum with a large enhancement at the unparticle threshold. 

The active-like mode features a stable isolated pole below the unparticle threshold, 
but "inherits" a non-trivial spectral density above it as a consequence of the mixing, 
even in absence of standard model interactions. The non- vanishing spectral density 
may open up new kinematic channels for weak interactions. 

• We obtain the radiative decay width of the unsterile-like mode into the active-like 
mode and a photon via a charged current loop. For large anomalous dimension (but 
7y < 1/3) we find a substantial suppression of the radiative decay width suggesting a 
concomitant weakening of the bounds on the mass and mixing angle from the X-ray 
or soft gamma ray background. 

II. UNSTERILE- ACTIVE MIXING 

In order to study the fundamental aspects of unsterile-active neutrino mixing and the 
potential cosmological consequences, we consider the simplest case with one unsterile and 
one active Dirac neutrino. The case of Majorana neutrinos and a triplet of unsterile neutrinos 
as envisaged in extensions beyond the standard model [10\ will be studied in detail elsewhere. 

^ Ref. [61] speculated that fermionic unparticle coupled to an active neutrino might give rise to energy- 
dependent mixing. 
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For an unsterile Dirac fermion the Lagrangian density in momentum space is [29l [59] 



C = iJu{-p){f^-M)F{p)iJu{p), 



(2.1) 



where 



F{p) 



A2 . 



; < ?7 < 1 . 



(2.2) 



A is the scale below which the low energy dynamics is dominated by the infrared fixed 
point of the conformal sector. Below this scale the unparticle is described by an interpo- 
lating field whose two point correlation function scales with an anomalous dimension [24j. 
Consistency of the unparticle interpretation requires that 



M < A. 



(2.3) 



We consider the mixing with an "active" massless Dirac neutrino ijja of the form 



(2.4) 



A see-saw mechanism consistent with the unparticle nature of the sterile neutrino, namely 
an interpolating effective field below a scale A entails the following hierarchy of scales 



m <C M < A. 



(2.5) 



In what follows we will explicitly invoke this hierarchy in the analysis. 
It is convenient to introduce the "flavor doublet" 



and write the Lagrangian density for the unsterile and active fermions as 



C = *(-j9) - M ^'(p) 



where 



and 



F = 



M 



1 

F{p) 

m 
m MF{p) 



(2.6) 



(2.7) 



(2i 



(2.9) 
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The equation of motion is 

- M = 0. 
It is convenient to introduce the spinor $ so that 

qf(p) = ^I + F-^M $(j9), 

which obeys the following equation of motion 

- M F"^ ]m1$(p) = 0, 



where 



M F"^ M = I ^^P^ 

mM + M'^F{p) 



The matrix in (2.12) can be written as follows 



/F - M F"^ M = a{p) I - p{p) 



-C{p) Sip) 
S{p) C{p) 



Introducing the shorthand 



we obtain 



Q = p'- 



m 

F{p} 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



and 



a{p) = 

m = 



Q + F{p){Q-M') 



Q-Fip){Q-M') 



C{p) 
S{p) 



Q-F{p){Q-M^ 



2 Pip) 



mM 



These functions obey 



C\p)+S\p) = l. 



(2.16) 
(2.17) 



(2.18) 
(2.19) 

(2.20) 



It becomes clear that the dispersions relations for the propagating modes correspond to 



det 



/F - M F"^ M = a^(p) - 0^ip) = ^ QiQ - M^) - ——— = 0, (2.21) 
J F[p) 



which then leads to a self-consistent equation for the dispersion relation of the propagating 
modes 



pI = 



1± Wl 



4m2 



(2.22) 



M2F(p±) 

The Klein-Gordon operator ( |2.14[ ), which we obtained from squaring the Dirac operator, 
can be diagonahzed using the transformation 



U{p) 



c{p) s{p) 
-s{p) c{p) 



where 



c{p) = ^[l+C{p)f ; s{p) = ^[l - C{p) 



The resulting spinor is given by 



(2.23) 



(2.24) 



(2.25) 



and in this basis, the Klein-Gordon matrix (2.14) becomes 



U-\p) p^F - M F"^ M Uip) 



(2.26) 



a{p)+p{p) 

aip) - p{p) 

An alternative manner to understand the solutions to the equations of motion is by going 
to the chiral representation. We do this by expressing the spinor ^' in terms of its right and 
left components, each a flavor doublet 



R 



(2.27) 



We can then expand the right and left components in the helicity basis 
where 



h=±i 



a • p 



v''{p} = h v''{p) ; h = ±l, 



and are flavor doublets. We find that 

(p° - %1)F + M = 0, 

(/ + %1)Fe^ + M^^ = 0. 



(2.28) 
(2.29) 

(2.30) 
(2.31) 



Using (2.31), we can expres ^'l in terms of and obtain 



F-1 M 



(po + h\p\) 



(2.32) 



Alternatively, we can express in terms of using (2.30), and obtain 

F-i M 



p^F - M F"^ M 
Introducing the "mass eigenstates" 



(2.33) 



(2.34) 



it follows that 



a{p) + I3{p) 






a{p) - Pip) 



U-\p)(f-'^ M)uip) 



ip° + h\p\) 



X^R- 



(2.35) 



or 



a(p)-/?(p)r' ' """^ 



(2.36) 



It is clear that although the transformation U diagonalizes the Klein-Gordon operator. 



for F 7^ 1, it does not diagonalize the flavor matrix F in (2.35) and (2.36). 

The dispersion relations of the propagating eigenstates correspond to the solutions of 



a{p) = ±P{p), which are determined by the self-consistent equation (2.22). The roots p± 
correspond to a{p±) = ±P{p±) respectively. 

The propagator for the flavor doublet ^, denoted by S, obeys 



^F-M %{p) = I, 



(2.37) 



where I is the identity in both flavor and Dirac space. Pre-multiplying (2.37) by + 
MF-i^ 



we obtain 



/F - M F"^ M S(p) = ^I + MF"^ 
In the new basis '^(p) = U~^{p) ^(p), the propagator is expressed by 



(2.38) 



U-\p)S{p)U{p) 









1 



^l + U-\p)[MF-^]U{p) 



(2.39) 



a(p)-l3(p) 
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Only for F = 1 is the matrix inside the bracket on the right hand side of (2.39) diagonal. For 



F 7^ 1, there is no unitary transformation that diagonalizes both the matrices proportional 
to ^ and MF-^ 

// F{p) is real, C{p),S{p) introduced in ( |2.14|2.18|2.l'9l ) can be identified as cosine and 
sine of (twice) of the mixing angle, namely 



c{p) = cos{2em{p)) ; s{p) = sm{2eup)), 

and similarly c{p), s{p) in ( 3.8|2.24 ) 

c{p) = cos(6'^(p)) ; s{p) = sm{Om{p)). 



(2.40) 



(2.41) 



Therefore, if F{p) is real, the transformation U{p) is unitary and 9m{p) is identified as 
the mixing angle. However, F{p) becomes complex above threshold p'^ > refiecting the 
multiparticle nature of the unparticle interpolating field i/ju- Thus, C,5,c and s cannot be 
interpreted as cosine and sine of (twice) the mixing angle. 

For F{p) = 1, the unparticle field is just an ordinary Dirac spinor field with canonical 
kinetic term. In this case, C, 5, c and s become independent of p, and they are given by 

M ^ . . 2m 



C = cos{29o) 



VM2 + 4m2 



; 5 = sin(26^o 



VM2 + 4m2 



and 



c = cos(^o) ; s = sin(^o)- 
Here, the angle Oq is the usual mixing angle for the see-saw mass matrix. 



(2.42) 



(2.43) 



Although the transformation U (3.8) diagonalizes the Klein-Gordon operator in the equa- 



tions of motion (2.26), it does not diagonalize the propagator or the Lagrangian in terms 



of the "mass eigenstates" . Furthermore, the solutions of the equations of motion in the 
transformed basis, namely the spinor x 

( [235p36l ), still has mixing; between them. This is 
because U~^{p)(f~^ M^U{p) is not diagonal. It follows that there is no unitary transfor- 
mation that diagonalizes the propagator. The problem in the diagonalization can be traced 
to the spinor nature of the unparticle field, where there are two independent structures in 
the effective action, the mass matrix term and the kinetic term ^ multiplied by F. There 
simply is no unitary transformation that diagonalizes simultaneously both the mass term 
M and the kinetic term F. An alternative explanation using Lorentz invariance argument 
is presented in the Appendix lAl 
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III. NON-UNITARY TRANSFORMATION: FIELD REDEFINITION 



As pointed out in the previous section because of the spinor nature of the field and 
the fact that the matrix coefficients of the kinetic term ^ and the mass matrix do not 
commute, there is no unitary transformation that diagonalizes the full propagator, even for 



real F{p). However, the equation for the propagator (2.37) suggests that the following set 
of transformations will lead to a diagonalization of the propagator. Let us introduce 



(3.1) 



By multiplying the equation (2.37) on the right by \/¥ and on the left by l/ylF one finds 
the following equation for S 



where 



M = ^ M ^ 



The mass matrix M can be written as 



/F{P) 

M 



M ^ M 
M = — 1 + — 
2 2 



1 + 



4m^ 



where 



Cip) = 
S{p) = 



1 + 



M^F{p) 



-C{p) S{p) 
S{p) C{p) 



M^F{p) 
2m 



M ,/F{p) 



1 + 



4m 



M^F{p) 



When F{p) is real 



C{p) = cos{2ip{p)) ; S{p)sm{2ip{p)), 

where f{p) is a mixing angle that depends on p"^. 

It is clear that now the propagator can be diagonalized by the matrix 



U{p) 



c{p) s{p) 
-s{p) c{p) 



(3.2) 



(3.3) 



(3.4) 



(3.5) 
(3.6) 

(3.7) 



(3.^ 



11 



where 



1 r ~ 1^ ~ 1 



1 - Cip) 



(3.9) 



The mass matrix is now diagonal 



U-\p) M U{p) = Md 



Mi{p) 
M2{p) 



(3.10) 



with 



Mi(p) 
M2(p) 

If -F(p) is real, it follows that 



M 
T 

M 



1- 1 + 



4m^ 



1+ 1 



M2F(p) / 
M^F{p) j 



(3.11) 
(3.12) 



c{p) = cos(^(p)) ; s{p) = sin(^(p)) 



The transformed propagator 



is given by 



%mU-\p) S Uip) 



l/+Ml{p) r. 

p^-Mf(p) ^ 

A iy+M2(p) 



(3.13) 



(3.14) 



(3.15) 



The transformation (3.1 ) has a natural interpretation in terms of a field redefinition. This 



can be inferred from the form of the kinetic term for the unparticle field ipjj f) F{p) ipu^ 
which suggests that F{p) can be interpreted as a momentum dependent wave function 
renormalization. Let us define the rescaled field as 



which along with 



forms the fiavor doublet 



V = 



(3.16) 



(3.17) 



(3.18) 



12 



With this field redefinition, the Lagrangian density becomes 



(3.19) 



No physics has been lost with this field redefinition as the correlation functions of the 
original unparticle field ifju may be obtained as follows. Let us introduce Grassman sources 
Ju coupled to the unparticle field in the Lagrangian density, namely 



(3.20) 



Upon field redefinition (3.16), the source terms become vuJj^F{p\ etc. Furthermore, at 



the level of the path integral, the field redefinition multiplies the measure by an overall field 
independent constant which cancels in all correlation functions. 

The necessity of a field redefinition to rescale to unity the coefficient of has been also 
recognized in Ref. [62j within the context of radiative corrections in the quark sector of the 
standard model. 

The full Lagrangian density can now be diagonalized by introducing the "mass basis" 
z/i, V2 as 



U{p) 



^1 



^2 



(3.21) 



where U{p) is given by Eq. (3.8). 

The dispersion relations are obtained from the respective Dirac equations for the mass 
eigenstates, 

[fl-M,]u, = ; z = l,2 (3.22) 
leading to the self-consistent equations 



m 



M2 



1 + 



1+1 + 



Am? 



(3.23) 



(3.24) 



These dispersion relations are exactly the same as Eq. (2.22) with pi^2 = P-,+ respectively. 



The reason for this is that the determinant (2.21), which determines the dispersion relation, 
is simply rescaled, namely 



det 



/F - M F"^ M 



= det 



= F{p) det 



(3.25) 
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The original unsterile and active fields are related to the mass eigenstates as 



(3.26) 
(3.27) 



Therefore, arbitrary correlation functions of the unsterile field ipu can be obtained from the 
propagators of the mass eigenstates z^i,z^2- Furthermore, since the unsterile field does not 
couple to any other field of the standard model, only the unsterile-like mass eigenstate 1^2 



can participate in weak interaction processes via the relation (3.26) and this field does not 



directly involve the field redefinition (3.16). 

For < M^, F{p) is real and the transformation U{p) is unitary. U{p) is determined 
by the mixing angles defined by ( 3.5[3.6|3.7 ). However, above the unparticle threshold 
p^ > M^, F{p) is complex and the matrix U{p) is not unitary. This is a consequence of the 
coupling to a continuum of states. A similar situation emerges in the theory of neutral meson 
mixing, where the absorptive part of the Wigner-Weisskopf Hamiltonian, which describes the 
quantum mechanics of neutral meson mixing, prevents a diagonalization of the Hamiltonian 
via a unitary transformation. This situation has been analyzed in detail in Refs. [63l 
[64] . In particular, Ref. [63j discusses the reciprocal basis that corresponds to fields that 
are transformed by a non-unitary transformation. A similar discussion appropriate to the 
quark sector of the standard model is given in Ref. [62j and a non-unitary transformation 
concerning time-reversal violation in the neutral kaon system can be found in Ref. [65j . The 
reader is referred to these references for a detailed discussion of the reciprocal (or dual or 
in-out) basis within the theory of neutral meson mixing. 

The same analysis in terms of the reciprocal (or dual) basis applies to the case under 
consideration for p^ > when F{p) becomes complex. 



A. Complex Poles and Spectral Densities for the Active-like Mode 

The dispersion relations of the propagating modes, are obtained from the complex poles 
of the propagator corresponding to pi 2 = ^1,2- Self-consistent solutions of the equations 
(|3j23||3^ are in general difficult to obtain analytically, however progress can be made in 



the relevant case m <C M and assuming self-consistently that rri^ /M^F{pi^2) <^ 1- 
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With this approximation we find 



m 



M4 F2(p) 



+ 



(3.28) 



therefore the self-consistent equation (3.23) for p\ becomes 

,4 



m 



M4 F2(pi) 

Anticipating self-consistently that -C we write 

M2 



+ 



(3.29) 



Fipi) 



A2 



-v r 



1 



M2 



(3.30) 



leading, to lowest order in the ratio m /M F, to the solution 



Pi 



m 
M2 



M2 

A2 



2»? 



(3.31) 



namely an isolated pole below the multiparticle threshold at p2 = j^f 2 Near this pole we 
find 

1 ^1 



where 



p^-M^{p) p^-Mf' 
M2 



^ 1 + 2?7 



M2 ■ 



(3.32) 



(3.33) 



For F{p) = 1 {i] = 0), Ml = m? /M is recognized as the smallest eigenvalue of the see-saw 
mass matrix, namely the mass of the lightest neutrino. This pole lies on the real p^ axis and 
describes a stable active- like propagating mode. 

The active-like propagator also features a discontinuity across the real axis in the complex 
p2-plane for p^ > M"^, since 



m 
M2 



-p2 + m2 - ie 
A2 



2»? 



(3.34) 



It is convenient to introduce the dimensionless variables 



X 



p^-M^ 
M2 



A = 2 



m 
M2 



M2 



A2 



(3.35) 



and use these to define the dimensionless spectral density 

p,{x) = — Disc 



(3.36) 
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where the discontinuity is non- vanishing for x > 0. We find that 



e{x) f x^^ sm{27rr]) 



7T 



1 2 



+ 



1 2 • 



sin(27r7y) 



(3.37) 



X + 1 — ^ x^^ cos(27r7y) 

This spectral density vanishes at threshold p'^ = (x = 0), increases rapidly reaching a 
broad maximum and diminishes for increasing x (Fig. ([T])). 



0.0020 



1^ 

a 



0.0015 - 



0.0010 



0.0005 



0.0000 




FIG. 1: Spectral density for the active-like mode. 



It is remarkable that in the absence of other interactions, the propagator of the active-like 
(lightest) mass eigenstate features a non-vanishing spectral density away from its mass shell 
for > > . This mode "inherits" a coupling to the continuum "hidden" sector as 
a consequence of the mixing with the unparticle. The non-vanishing spectral density above 
the unparticle threshold at p'^ = may lead to opening new kinematic channels when the 
active-like neutrino is coupled to the standard model fields. 



B. Complex Poles and Spectral Densities for the Unsterile-like Mode 



For m <C M, the self-consistent equation (|3.24|) for the unsterile-like mode becomes 

(3.38) 



-pI + -ie 



A2 



+ 



In terms of the dimensionless variables (3.35) this equation becomes 

1 V 



x = A 



X — 16 



(3.39) 
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We find that there is a solution only for 

Re{x) > , 0<r] < 1/3, 

and it is given by 



pI = 1 + Ai-^ I cos - ^sin 



TTT] 

I sm I 

rjy \1 — T] 



(3.40) 



This solution describes a pole in the complex plane (a resonance) and near this pole we find 



1 



where 



p2 - M|(p) - M| + iM2T ' 

TTT] 



Ml 



1 + Ai-'' cos 



and^ 



M2 1 

Ai-') sm , 

Ma VI 



1 



1 

7rr] 



?7 



7] 



1—7] 



(3.41) 
(3.42) 
(3.43) 
(3.44) 



The imaginary part is a consequence of the fact that the real part of the pole is above the 
unparticle continuum determined by the multiparticle threshold at = M^. For F{p) = 1 
(r; = 0), the largest eigenvalue of the see-saw mass matrix is at M + m^/M > M. After 
mixing, the new pole is in the unparticle continuum, moving off the real axis into a second 
(or higher) Riemann sheet in the complex p'^ plane. The imaginary part describes the decay 
of the unsterile like mode into the active-like mode and particles in the "hidden" conformal 
sector. A similar phenomenon was observed in the bosonic case in Ref. [48]. We refer to F as 
the "invisible width" of the unsterile-like neutrino since it describes its decay into an active- 
like and conformal massless particles in the "hidden sector." A more detailed discussion and 
interpretation of this result, based on a renormalization group |S2l EH EH] resummation is 
presented in Sec. |III C 



The spectral density is obtained from the discontinuity across the real axis in the complex 

A sin(7r7y) 



p'^ plane (see Eq.(3.36)) 



P2{x) 



e(x) 



X — Ax^ cos(7r7y) 



1 2 



Ax^ sin(7r7y) 



1 2 ' 



(3.45) 



^ We note that the unparticle field is not canonical, therefore the residue at the unparticle-like "pole" is 
not restricted by canonical commutation relations to obey < Z2 < 1. 
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which is displayed in Fig. ([2]). For < ?7 < 1/3, there is a resonance with a maximum 
confirmed to be given by 

(3.46) 



Xp = cos 



as obtained in Eq. (3.42) 



A = 0.3 



1 — 7] 




r| = 0.1, Xp = 0.247 

Ti = 0.15, Xp = 0.206 
T| = 0.20, Xp = 0.157 
= 0.25, Xp = 0.1 



FIG. 2: p2(x) vs. x for A = 0.3 and several values of 77. Xp is the real part of the complex pole 
( [3l2pl6l ). 



For A <C 1 the real part of the unsterile-like pole is very near the threshold at = 
and therefore the dimensionless ratio 



2tan(f^ 

' 1— ?7 



M2-M 



1 + A^ cos 



(3.47) 



determines if the resonance is broad or narrow as compared to the distance between threshold 
and the resonance. 

For Tj {F{p) —> 1), the width becomes very small and the resonance is sharp, centered 
at the mass of the sterile neutrino M + jM. As 77 ^ 1/3 from below, the real part of 
the pole approaches threshold [x^ 0) while the width remains constant. The resonance 



broadens enormously since the ratio (3.47) diverges as 77 ^ 1/3 , and merges with the 



threshold at = 1/3. There are no solutions of the self-consistency condition (3.38) for a 
complex pole for 77 > 1/3. 
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C. Resummation Interpretation of the Decay Width 



Consider the unparticle to be a single Dirac fermion (neutrino) interacting with massless 
particles of the conformal hidden sector, described by a conformal field A. Let us assume 
such interaction to be of the form 



Ant = ^*^*, 



(3.48) 



where 5^ is a small dimensionless coupling. In perturbation theory the self-energy of the 
''unparticle field" is depicted in Fig. ^ 




11/ 



FIG. 3: ^ self-energy to lowest order in g. A is the conformal field. 



To lowest order in g, the self energy, which is once-subtracted to vanish at ^ = M, is 
given near the mass shell by 

p"^ + M'^ -ie 



E(p) = -r]{^ - M) In 
where A is a renormalization scale, 



A2 



(3.49) 



?7 = cg\ (3.50) 

and c is a constant that depends on the nature of the conformal field A (gauge or scalar 
massless particle). 

Integrating out the conformal field A leads to the following effective action for ^ 



*(^- M) 



1 — ?7 In 



A2 



+ - ie 
A2 



(3.51) 



where in the last line we have invoked a renormalization group resummation [521 [541 [55] of 
the infrared threshold divergences. The infrared logarithmic divergence at = and the 
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imaginary part for > is the result of the emission and absorption of massless quanta, 
an ubiquitous phenomenon in gauge theories (for a discussion within QCD, see [52]). 

Now consider couphng the heavy neutrino to a massless (active) neutrino ijja via the 
coupling 

Crri = ^ rnipa + h.c. , (3.52) 



leading to a see-saw mass matrix of the same form as in Eq. (2.9) for F{p) = 1 with m <C M. 
The see-saw mass matrix can be diagonalized by a unitary transformation with the mixing 



angle Oq^ given by the relations (2.42). The fields that describe the "mass eigenstates" are 



V^i = cos(6>o)'0a + sin((9o)^^, 
7/^2 = cos((9o) * - sin((9o) ipa , 

where the masses corresponding to the fields ^01^2 are Mi^2 respectively, with 

2 2 



To lowest order in the see-saw ratio m/M, it follows from (2.42) that 

m 



sin(6'o) 



M 



(3.53) 
(3.54) 

(3.55) 

(3.56) 



The — ipa mixing leads to the following interaction vertex between the fields associated 
with the mass eigenstates and the conformal field A 

(3.57) 



^int = g (^cos((9o)'02 - sin((9o)'0ij^(^cos(6^o)'02 - sin((9o)'0i 
The self-energy for the 7/^2 field now includes the diagram depicted in Fig.([4]). 




FIG. 4: Self-energy for mass eigenstate '02- 



The cut discontinuity across the intermediate state relates the absorptive (imaginary) 
part of the self-energy on the mass shell of the external fermion to its decay rate. This 
relation is depicted in Fig. ([s]). 
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FIG. 5: Cutkosky cut for the self-energy for mass eigenstate ^2- The imaginary part of the self- 
energy determines the width F in the pole of the propagator. 



A standard, straightforward calculation of the decay rate for the process ip2 
taking ipi to be a massless scalar and Dirac fermion, respectively, yields 



27rr] M2 ( sin(^o) cos(0o) ) ^ 



2m^ 



(3.58) 



where 77 is given by Eq. (3.50) and we have used the approximations M2 ^ M; sin(^o) 



m/M to lowest order in the see-saw ratio m/M <C 1. 

This result coincides to lowest order in r] and m/M with the non-perturbative imaginary 
part of the unparticle-like pole given by Eq. ( 3. 43] ). 

This simple analysis confirms that the imaginary part F in the propagator of the 



unparticle-like mode (3.41) describes the decay of the unparticle-like mode into the active- 



like mode and particles in the hidden conformal sector. This analysis also validates the 



interpretation of the width of the unsterile mode (resonance) F given by (3.43) as an "in- 



visible width" as opposed to the radiative decay width that arises via weak interactions 
described below. 



IV. COSMOLOGICAL CONSEQUENCES: RADIATIVE DECAY OF THE 
UNSTERILE-LIKE NEUTRINO AND THE X-RAY BACKGROUND 

Although sterile neutrinos only couple to active neutrinos via an oS diagonal mass matrix, 
the diagonalization of this mass matrix results in effective couplings between the sterile-like 
neutrino mass eigenstate and standard model particles, namely active neutrinos and charged 
leptons. Consider the simple case of (canonical) sterile neutrinos coupled to active neutrinos 



via a see-saw mass matrix of the form (2.9) with F = 1, diagonalized by the usual unitary 
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transformation. In this case, 



= cos(0o) ^1 + sin(^o) ^2 : 
= cos(^o) ^2 — sin(0o) ^1 



(4.1) 
(4.2) 



where as usual z^i, V2 are the hght (active-hke) and heavy (sterile-hke) neutrino mass eigen- 
states, with masses Mi ?^ w? jM\ M2 ^ M -\- /M ^ Mi, respectively. The mixing angle 
is determined by Eqs. ( 3.7[3.5|3.6 ) with F = 1. 

The charged current interaction yields an interaction between the sterile-like neutrino 
and the charged lepton 



Ccc = g^aL )^/L = ^(cos((9o)z^iL + sin(6>o)z^2L) Wh- 



(4.3) 



This interaction vertex leads to the radiative decay of the sterile-like neutrino z/2 ^ ^^i 7 
The diagrams that describe this process in unitary gauge are shown in Fig. ([6]). For 
M2 ^ Ml, the radiative decay width is given by [15] 



r ~ 



32^ 



mi 



1 4 



w 



sm^{9o) cos2(^o) 



(4.4) 



where mi is the mass of the charged lepton in the loop in Fig. ([g]). 




W 




^2 / 



^1 



FIG. 6: Contributions to radiative decay of a sterile-like neutrino in unitary gauge. 



If sterile neutrinos are suitable dark matter candidates, with M2 of order of a few keV 
[3 El O [ini [m [ini [H] , the radiative decay of the sterile-like neutrino yields a contribution 
to the X-ray or soft gamma ray background, from which stringent bounds on the mass and 
mixing angle of the sterile-like neutrino are obtained [61 [lOl 1121 113 1^ [19] . 

The calculation leading to the radiative decay width ( |4.4 ) uses the standard assumptions, 
namely that the propagator of 1^2 is that of a free particle featuring a single particle pole. 
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that the mixing angles are independent of momentum, and that the in and out states are 
in their (single particle) mass shells. All of these assumptions must be revised in view of 
the results obtained in the previous sections. We will proceed to obtain an estimate of the 
unparticle effects upon the radiative decay width by modifying the calculation leading to 



(4.4) by including the following unparticle effects: 



We will consider the mass eigenstate V2 to be described by the propagator (3.41) but 
neglecting the "invisible width" F. This is similar to the situation in calculating the 
decay of a vector boson by considering it to be an asymptotic state with a propagator 
featuring a single particle pole. This assumption restricts the validity of our estimate 
to0<7y<l/3 since for 77 > 1/3, the spectral density of 1^2 is a broad continuum 
above threshold. Furthermore, since the residue at the pole ^ M| is a finite 
wave function renormalization Z2, the total transition probability will be multiplied 
by Z2 = 1/(1 — ry). Similarly, the mass eigenstate Vi is described by the propagator 
(3.32), which does feature an isolated single particle pole at = Ml but with residue 



Zi = 1/(1 + r]A'^/2y Therefore, the transition probability must also be multiplied by 
the wave function renormalization Zi. 



The mixing angles cos(0o) ^(p); sin(^o) s{p) are given by Eqs. ( 3.9[3.5[3.6 ), which 



depend on the momentum. The radiative decay rate corresponds to setting both the 
decaying and the product particles on their mass shells. Therefore, c{p)^ which is the 
amplitude of Ui in must be evaluated on the mass shell of the active-like mode, 
namely = , while ^{p)^ which is the amplitude of U2 in z^^, must be evaluated at 
the mass shell of the unsterile-like mode, p'^ = M|. Because is below threshold, it 
follows that c{p) is real, however, since M| is above threshold, ^{p) is complex. Being 
a probability, the decay rate involves the modulus squared of these quantities, namely 



in the expression (4.4), we must replace 

1 



2 



cos^(^o) ^ c\p) = - 1+C{p) ^ ^, (4.5) 



sin^eo) ^ \^ip)\ = ll-C{p) ^ ^. (4.6) 



2 



For 1] = {F{p) = 1) and m <C M, it follows that cos(^o) ~ 1, sin(^o) ~ m/M. 
In the same limit, we find that for F{p) 7^ 1, c{p^ = Ml) f» 1, s(p^ = M|) fa 
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m 



/M\jF{p^ = M|) (see Eqs. ( 3.5|3.6 )). The overall change then corresponds to 

1 



cos'(^o)sin^(^o) 



F{p) 



cos^(^o)sin^(^o), 



(4.7) 



where p'^ = pi is given by (3.40) 



Including all of these modifications we obtain the ratio of the radiative decays for the 
unsterile, and the (canonical) sterile neutrino as 



A ^ 



A2 



r.™ {i-v){i + vf)' 

For a see-saw mass matrix with m <^ M, we recognize that 

A fti 2sin^(6'o) 



A2 



(4.8) 



(4.9) 



where is the mixing angle for r] = (namely the case of canonical sterile neutrinos mixing 



with active ones). Therefore, we can write the ratio (4.8) as 



1-77 



(4.10) 



r.,^.,, (1 - 77) (1 + vf) 

Taking sin^ <C 1 and M < A, consistent with a large see-saw and with the unparticle 
interpretation of the sterile neutrino below a scale A, respectively, we see that the unparticle 
nature of the sterile neutrino can lead to a substantial suppression of the radiative decay 
rate. Even for 77 <C 1, the current bounds on the mixing angles and masses for sterile 
neutrinos from the observations of the X-ray or soft gamma ray background can be weakened 
considerably. As an example, taking sin(0o) ^ 10~^, M ~ keV and A ~ TeV, which are 
within the range of expectation for physics beyond the standard model, and taking ~ 0.1 



as an example inspired by results in QCD [52j^, we find that the ratio (4.10) < O{10~^). 

This is one of the main cosmological consequences of the unparticle nature of the sterile 
neutrino. 



^ The value of rj for QCD obtained in Ref. [52], 77 ~ 0.5 is larger than 1/3 which hmits the validity of our 
result assuming a sharp resonance for We have taken 7^ ~ 0.1 as a representative value of the range 
for large anomalous dimensions in QCD. 
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V. CONCLUSIONS AND FURTHER QUESTIONS 



In this article we considered the possibihty that the SU (2) singlet sterile neutrino might 
be an unparticle, an interpolating field that describes a multiparticle continuum as a conse- 
quence of coupling to a "hidden" conformal sector and whose correlation functions feature 
an anomalous scaling dimension 77. We studied the consequences of its mixing with an active 
neutrino via a see-saw mass matrix. We focused on the simplest setting of one unsterile and 
one active Dirac neutrino, postponing a more detailed study of Majorana neutrinos and 
several flavors to further study. Our goals here are two-fold: 

1. to study the consequences of the mixing between the non-canonical unsterile with 
a canonical active neutrino, along with the corresponding dispersion relations and 
propagating states, 

2. to explore cosmological consequences, in particular the radiative decay width of the 
unsterile-like neutrino into the active-like and a photon, via charged current loops, and 
to establish how the unparticle nature of the sterile neutrino modifles the radiative 
line-width, which is an important tool to constrain the mass and mixing angles of 
cosmologically relevant sterile neutrinos. 

We found that the mixing between a non-canonical and a canonical fermion fleld exhibits 
several unexpected subtleties. There is no unitary transformation that diagonalizes the full 
propagator due to the non-canonical nature of the unsterile neutrino. This forces us to make 
a field redefinition for its complete diagonalization. 

The unsterile-like propagating mode is described by a complex pole above the unparticle 
threshold for < 77 < 1/3, featuring an "invisible width". A perturbative analysis and a 
renormalization group inspired resummation suggest that this width results from the decay 
of the unsterile-like mode into an active-like mode and particles in the conformal sector. As 
77 l/3~, the complex pole merges with the unparticle threshold and disappears, while 
the spectral density of the unsterile neutrino features a broad continuum above threshold 
with a threshold enhancement. The active-like mode corresponds to a stable particle, whose 
propagator features an isolated real pole below the unparticle threshold. This mode "inher- 
its" a non-vanishing spectral density above this threshold, even in the absence of standard 
model interactions. This novel feature may potentially have relevant consequences since the 
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non-vanishing spectral density may open new kinematic channel for standard model pro- 
cesses even to lowest order in weak interactions, a possibility that will be studied in detail 
elsewhere. 

Considering unsterile neutrino as a dark matter candidate, we studied the influence of 
the unparticle nature of the sterile neutrino on the radiative decay width into an active 
neutrino and a photon via charged current loops. We flnd the ratio of decay widths between 
the unparticle case and the canonical case to be 



pf7 



2sin2(^o 



' A2 



1-77 

(5.1) 



r.,^.,^ (1 - 77) (1 + rif) ' 

where is the mixing angle for = 0, M is approximately the mass of the unsterile like 
neutrino and A ^ M is the unparticle scale. This ratio suggests a substantial suppression 
of the radiative decay line width for M ^ keV and A ^ TeV, even for < 0.1. This results 
in a weakening of the bounds on the mass and mixing angle from the X-ray and soft gamma 
ray backgrounds. 

Of course, a detailed assessment of the suppression of the radiative decay width hinges 
on the (unknown) values of A and ry, which may emerge from the experimental program at 
LHC in the exploration of physics beyond the standard model. 

To further explore the possibility of an unsterile neutrino as a dark matter candidate, 
understanding its production process is necessary. Since an unsterile neutrino only interacts 
directly with the active one, the most effective dark matter production mechanism in this 
scenario is via unsterile-active neutrino oscillations. It would be interesting to study the 



implications of our results in Section [TTT] in the dark matter production mechanism along 
the line of Ref . [2j. 
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APPENDIX A: AN ALTERNATIVE EXPLANATION FOR THE NON- 
UNITARY TRANSFORMATION 



In the helicity basis of Eq. (2.28), the Lagrangian density in momentum space is given 



by 



(/ - h\p\)¥ M 

M (pO + h\p\)F 




(Al) 



Lorentz invariance does not allow us to mix the (2, l)-representation of the Lorentz group 
with the (1, 2)-representation. Therefore, to diagonalize the action, the allowed transforma- 
tion is ^^ j^ ^R,lCr,l^ such that all the following matrices: 



U 



R 



{p^-h\p\)F Ul (/ + /i|p1)F Ul\ UrMUI^ and UlMU^^ (A2) 



are all diagonal. Since F is diagonal, and yet not proportional to the identity, the only 
possible unitary transformations that diagonalize the first two are 



URr = 



1 

±1 



or 



±1 
1 



(A3) 



However, none of the combinations of these possibilities diagonalize the last two matrices in 



(A2). Therefore, there is no unitary transformation that diagonalizes the full propagator. 



[1] For a thorough description of the see-saw mechanism and list of references see Massive neu- 
trinos in physics and astrophysics^ R. N. Mohapatra, P. B. Pal, (World Scientific, Singapore, 
1997). 

[2] S. Dodelson and L. M. Widrow, Phys. Rev. Lett. 72, 17 (1994) |arX iv:hep-ph/9303287| . 
[3] S. Colombi, S. Dodelson and L. M. Widrow, Astrophys. J. 458, 1 (1996) |arXiv:astro- 
|ph/9505029 . 

[4] X. D. Shi and G. M. Fuller, Phys. Rev. Lett. 82, 2832 (1999) | arXiv:astro-ph/9810076 |; 
K. Abazajian, G. M. Fuller and M. Patel, Phys. Rev. D 64, 023501 (2001) |arXiv:astro- 
|ph/0101524j ; K. N. Abazajian and G. M. Fuller, Phys. Rev. D 66, 023526 (2002) |arXiv: astro- 
|ph/0204293 |; G. M. Fuller, A. Kusenko, I. Mocioiu and S. Pascoh, Phys. Rev. D 68, 103002 

27 



(2003) farXiv: astro-ph /03072671 ; K. Abazajian, Phys. Rev. D 73, 063513 (2006) |arXiv:astro- 
|ph/05 12631]. 

[5] M. Shaposhnikov and I. Tkachev, Phys. Lett. B 639, 414 (2006) | arXiv:hep-ph/0604236| ; 
M. Shaposhnikov, JHEP 0808, 008 (2008) |a rXiv:0804.4542 [hep-ph]] , |arXiv:astro- 



|ph/0703673[ 

[6] A. Kusenko, AIP Conf. Proc. 917, 58 (2007) |ar Xiv:hep-ph/0703116] , Int. J. Mod. Phys. D 



16, 2325 (2007) | arXi v:astro-ph/0608096| ; T. Asaka, M. Shaposhnikov and A. Kusenko, Phys. 
Lett. B 638, 401 (2006) fa rXiv:hep-ph/0602150| ; P. L. Biermann and A. Kusenko, Phys. Rev. 



Lett. 96, 091301 (2006) |arXiv:astro-ph/0601004|. 



[7] K. Petraki and A. Kusenko, Phys. Rev. D 77, 065014 (2008) |arXiv:071 1.4646 [hep-ph 



[8] K. Petraki, Phys. Rev. D 77, 105004 (2008) |arXiv:0801.3470 [hep-ph 



[9] D. Boyanovsky, H. J. de Vega and N. Sanchez, Phys. Rev. D 77, 043518 (2008) 
|arXiv:0710.5180 [astro-ph ]]. 
[10] A. Boyarsky, J. W. den Herder, A. Neronov and O. Ruchayskiy, Astropart. Phys. 28, 



303 (2007) |arXiv:astro-ph/0612219|, A. Boyarsky, D. lakubovskyi, O. Ruchayskiy and 



V. Savchenko, Mon. Not. Roy. Astron. Soc. 387, 1361 (2008) |i arXiv:0709.2301 [astro-ph]| , 
A. Boyarsky, D. Malyshev, A. Neronov and O. Ruchayskiy, Mon. Not. Roy. Astron. Soc. 387, 



1345 (2008) [arXiv:0710.4922 [astro-ph]|, A. Boyarsky, O. Ruchayskiy and M. Shaposhnikov, 



|arXiv:0901.0011 [hep-ph] 



[11] D. Boyanovsky, Phys. Rev. D 78, 103505 (2008) |arXiv:0807.06 46 [astro- ph]]. 
[12] M. Loewenstein, A. Kusenko and P. L. Biermann, |arXiv:0812.2710 astro-ph ] 



[13] H. J. de Vega and N. G. Sanchez, arXiv:0901.0922 [astro-ph.CO 



[14] J. Wu, C. M. Ho and D. Boyanovsky, |arXiv:0902.4278 [hep-ph]] 



[15] P. B. Pal and L. Wolfenstein, Phys. Rev. D 25, 766 (1982). 

[16] A. D. Dolgov and S. H. Hansen, Astroparticle Physics 16, 339, (2002) |arXiv:hep-ph/0103118[ 



[17] S. Riemer-Sorensen, S. H. Hansen and K. Pedersen, Astrophys. J. 644, L33 (2006) 
[arXiv: astro- ph/06036 61 1 ; S. Riemer-Sorensen, K. Pedersen, S. H. Hansen and H. Dahle, Phys. 



Rev. D 76, 043524 (2007) [ar Xiv:astro-ph/0610034] . 
[18] . N. Abazajian, M. Markevitch, S. M. Koushiappas and R. C. Hickox, Phys. Rev. D 75, 063511 



(2007) [ar Xiv:astro-ph /06 1 1 144] . 
[19] S. Riemer-Sorensen and S. H. Hansen, arXiv:0901.2569 [astro-ph]] 



28 



[20] V. G. Sinitsyna, M. Masip, S. I. Nikolsky and V. Y. Sinitsyna, [arXiv: 0903 .4654 [hep-ph 



[21] A. A. Aguilar-Arevalo et al. [The MiniBooNE Collaboration], Phys. Rev. Lett. 98, 231801 
(2007) [arXiv:0704.1500 [hep-ex] |. 



[22] M. Maltoni and T. Schwetz, Phys. Rev. D 76, 093005 (2007) [arXiv:0705.0107 [hep-ph; 



[23] S. N. Gninenko, |arXiv:0902.3802 [hep-ph][ 



[24] H. Georgi, Phys. Rev. Lett. 98, 221601 (2007) |arXiv:hep-ph/ 0703260]; H. Georgi, Phys. Lett. 
B 650, 275 (2007) [arXiv:0704.2457 [hep-ph]]. 

[25] K. Cheung, W. Y. Keung and T. C. Yuan, Phys. Rev. Lett. 99, 051803 (2007) | arXiv:0704.2588. 
I [hep-ph] j ; A. Arhrib, K. Cheung, C. W. Chiang and T. C. Yuan, Phys. Rev. D 73, 075015 
(2006) [arXiv: hep-ph /0602175|; for a recent review see: K. Cheung, W. Y. Keung and 
T. C. Yuan, AIP Conf. Proc. 1078, 156 (2009) [arXiv:0809.0995 [hep-ph] |. 



[26] J. L. Feng, A. Rajaraman and H. Tu, Phys. Rev. D 77, 075007 (2008) [arXiv:0801.1534 [hep- 



|ph] |; A. Rajaraman, AIP Conf. Proc. 1078, 63 (2009) |a rXiv:0809.5092 [hep-ph]] . 
[27] A. Delgado, J. R. Espinosa and M. Quiros, JHEP 0710, 094 (2007) [arXiv:0707.4309 [hep- 



|ph] |; A. Delgado, J. R. Espinosa, J. M. No and M. Quiros, Phys. Rev. D 79, 055011 (2009) 
|arXiv:08r2.1170 [hep-ph]| . 
[28] T. Banks and A. Zaks, Nucl. Phys. B 196, 189 (1982). 



[29] P. J. Fox, A. Rajaraman and Y. Shirman, Phys. Rev. D 76, 075004 (2007) | arXiv:0705.3092 
[[hep-ph]]. 



[30] M. Luo and G. Zhu, Phys. Lett. B 659, 341 (2008) |arXi v:0704.3532 [hep-ph ] |; Y. Liao, Phys. 

Lett. B 665, 356 (2008) |arXiv:0804.0752 [hep-ph]| ; R. Basu, D. Choudhury and H. S. Mani, 

[arXiv: 0803.4110 [hep-ph], 
[31] B. Grinstein, K. A. Intriligator and L Z. Rothstein, Phys. Lett. B 662, 367 (2008) 



|ar Xiv:0801.1140 [hep-ph]| 



[32] Y. Liao and J. Y. Liu, Phys. Rev. Lett. 99, 191804 (2007) [arXiv:0706.1284 [hep-ph]] 



[33] S. Hannestad, G. Raffelt and Y. Y. Y. Wong, Phys. Rev. D 76, 121701 (2007) | arXiv:0708.1404 
I [hep-ph ]]. 



[34] N. G. Deshpande, S. D. H. Hsu and J. Jiang, Phys. Lett. B 659, 888 (2008) |arXiv:0708.2735 



hep-ph]j . 



[35] A. Freitas and D. Wyler, JHEP 0712, 033 (2007) [ar Xiv:0708.4339 [hep-ph]] . 



[36] J. McDonald, |arXiv:0805.1888 [hep-ph]] JCAP 0903, 019 (2009) [a rXiv:0709.2350 [hep-ph 

29 



[37] H. Davoudiasl, Phys. Rev. Lett. 99, 141301 (2007) | arXiv:0705.3636 [hep-ph]] . 
[38] I. Lewis, |arXiv:0710.4147 [hep-ph]l 

[39] G. L. Alberghi, A. Y. Kamenshchik, A. Tronconi, G. P. Vacca and G. Venturi, Phys. Lett. B 
662, 66 (2008) |arXiv: 0710.4275 [hep-th]] . 



[40] S. L. Chen, X. G. He, X. P. Hu and Y. Liao, Eur. Phys. J. C 60, 317 (2009) |ar Xiv:0710.5129 
[[hep-ph]| . 



[41] T. Kikuchi and N. Okada, Phys. Lett. B 665, 186 (2008) |arXiv:071 1.1506 [hep-ph]]. 



[42] Y. Gong and X. Chen, Eur. Phys. J. C 57, 785 (2008) [arXiv:0803.3223 [astro-ph 



[43] H. CoUins and R. Holman, Phys. Rev. D 78, 025023 (2008) |arXiv:0802.4416 [hep-ph]] . 
[44] D. C. Dai and D. Stojkovic, |a'rXiv:0812.3396 [gr-qc]| 



[45] R. Zwicky, J. Phys. Conf. Ser. 110, 072050 (2008) [arXiv:0710.4430 [hep-ph]] 



[46] Y. f. Wu and D. X. Zhang, arXiv:0712.3923 [hep-ph]; R. Mohanta and A. K. Giri, Phys. Rev. 



D 76, 057701 (2007) [arXiv:0707.3308 [hep-ph]]; J. K. Parry, [ajXi^:0810.0971 [hep-ph] 



[47] G. J. Ding and M. L. Yan, Phys. Rev. D 78, 075015 (2008) [arXiv:0706.0325 [hep-ph]]; G. Bhat 



tacharyya, D. Choudhury and D. K. Ghosh, Phys. Lett. B 655, 261 (2007) |ar Xiv:0708.2835 
|[hep-ph]| . 

[48] D. Boyanovsky, R. Holman and J. A. Hutasoit, Phys. Rev. D 79, 085018 (2009) 

|ar Xiv:0812.4723 [hep-ph] |. 
[49] M. A. Stephanov, Phys. Rev. D 76, 035008 (2007) |a rXiv:0705.3049 [hep -ph]]. 
[50] D. Amit, Field Theory, the Renormalization Group and Critical Phenomena (McGraw-Hill, 

N.Y. 1978). 

[51] E. Brezin, J. C. Le Guillou, J. Zinn- Justin, in Phase Transitions and Critical Phenomena, 
Vol. 6 (Ed. C. Domb, M. S. Green, Academic Press, London, 1976) (page 127). 



[52] M. Neubert, Phys. Lett. B 660, 592 (2008) [a rXiv:0708.0036 [hep-ph]] . 
[53] F. Bloch and A. Nordsieck, Phys. Rev. 52, 54 (1937). 

[54] N.N.Bogoliubov, D.V.Shirkov, Introduction to the theory of quantized fields, (Interscience Pub- 
lishers, N.Y. 1959). 



[55] D. Boyanovsky and H. J. de Vega, Annals Phys. 307, 335 (2003) [arXiv:hep-ph/0302055| . 
[56] X. Q. Li, Y. Liu and Z. T. Wei, Eur. Phys. J. C 56, 97 (2008) [arXiv:0707.2285 [hep-ph]]. 



[57] S. Zhou, Phys. Lett. B 659, 336 (2008) [a rXiv:0706.0302 [hep-ph]] . 

[58] G. Gonzalez-Sprinberg, R. Martinez and O. A. Sampayo, Phys. Rev. D 79, 053005 (2009) 



30 



[arXiv:0808.1747 [hep-ph]] . 



[59] G. Cacciapaglia, G. Marandella and J. Terning, JHEP 0902, 049 (2009) |arXiv:0804.0424 
i [hep-p h]]. 



[60] G. von Gersdorff and M. Quiros, larXiv:0901.0006 [hep-ph; 



[61] T. Schwetz, JHEP 0802, Oil (2008) [ arXiv:0710.2985 [hep-ph]] . 

[62] Q. Duret, B. Machet and M. I. Vysotsky, Mod. Phys. Lett. A 24, 273 (2009) |arXiv:0810.4449 

|[hep-ph]|; |arXiv:0805.4121 [hep-ph]l 
[63] J. P. Silva, Phys. Rev. D 62, 116008 (2000) | ar Xiv:hep-ph/0007075] ; G. C. Branco, L. Lavoura, 

J. P. Silva, CP Violation, (Oxford University Press, Oxford, U.K., 1999). 



[64] M. Beuthe, G. Lopez Castro and J. Pestieau, Int. J. Mod. Phys. A 13, 3587 (1998) [a rXiv:hep- 
|ph/9707369l . 

[65] L. Alvarez-Gaume, C. Kounnas, S. Lola and P. Pavlopoulos, Phys. Lett. B 458, 347 (1999) 
|arXiv:hep-ph/9812326] . 



31 



